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ABSTRACT 

The de Broglie-Bohm ontological interpretation of quantum theory |2[] 
clarifies the understanding of many otherwise counter-intuitive quantum 
mechanical phenomena. We report here on an application of Bohm's 
quantum potential to the bonding and reactivity of small molecules. In 
the field of quantum chemistry, Bader has shown |]^ that the topology of 
the Laplacian of the electronic charge density characterizes many features 
of molecular structure and reactivity. Examination of high accuracy ab- 
initio solutions for several small molecules suggests that the Laplacian of 
Bader and the quantum potential of Bohm are structurally equivalent. 
It is possible that Bohmian mechanics using the quantum potential can 
make quantum chemistry as clear as it makes non-relativistic quantum 
mechanics. 
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A set of ab-initio computational experiments were performed simulating H2O, 
H2O2, and C2H4 to obtain the Laplacian of the electronic charge density (V^p). 
Bader has shown Q that V^p localizes lone pairs, bonded charge concentrations and 
regions subject to electrophilic or nucleophilic attack. V^p acts as an objective 
electron localization function^, ^, ^, |^ , and is free from the difficulties encountered 
when analyzing a wavefunction by decomposition into a particular orbital basis. 

In order to explore the application of the de Broglie-Bohm interpretation of 
quantum mechanics 0| to molecular physics, we also calculated the Bohm quantum 
potential Q for the same molecules. Visual comparison of the three dimensional 
structure of these two fields {Q and V^p) strongly suggests they are equivalent. 

Density functional theory calculations were carried out for the optimized ground 
state equilibrium geometries of water, hydrogen peroxide, and ethylene using Gaus- 
sian 94[|[ at the B3LYP/6-311++G(2d,2p) level of theory^. Isosurfaces of the Lapla- 
cian of the electronic charge density, V^p = C, were visualized [|T^] for these molecules 
over a range of values of C and are shown in figures 1-4 for the case of water [|TI | . 

Following the analysis of Bader P], regions of V^p < correspond to electronic 
charge concentration, and conversely regions of V^p > correspond to regions of 
charge depletion. The four views of H2O shown in Figure 1 correspond to a very 
negative value of V^p. The electronic charge concentrations associated with the 
nuclei as well as the two lone pair charge concentrations on the oxygen atom are 
visible. For a less negative value of the Laplacian (figure 2), concentrations of 
electronic charge on the OH bonds become visible. At V^p = (figure 3) the 
outline of the entire molecule emerges. Finally, for a positive value of V^p (figure 4), 
three regions of valence charge depletion are visible, as well as the region of depletion 
surrounding the Oxygen's core electrons. 

The structure of V^p may be formally defined in terms of the attractors, repellors, 
and saddle points in VV^p and their interconnections 0. In the case of water, there 
is one attractor at each nucleus, one in each lone pair region, one on each bond, and 
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no others (figures 1-3). 

To compute Bohm's quantum potential[|^ we first express the wavefunction in 
polar form, \1/ = Re'-^/^\ The quantum potential for an electron molecule using 
the Born-Oppenheimer approximation is then: 

Q is a function of the 3A^ coordinates of all of the electrons, and VQ is, in Bohm's in- 
terpretation of quantum mechanics, the nonlocal quantum force that acts in addition 
to the classical forces on each electron. 

For a one-electron molecule, Q reduces to 

ft ^ 

Visualization in three-dimensional physical space of the 3A^ dimensional quantum 
potential for multiple electron molecules was accomplished by calculating Q as 



a function of = yp/A^ where \E'i is the solution of the exact one-electron 
Schrodinger equation for an A^-electron molecule|jl2|. In this case, Q is identical 



to Qi times a constant. In terms of the charge density (p) and the Bader Laplacian 

g = vV2p-(Vp/2p)2. 

Isosurfaces of the (one electron) Bohm quantum potential for H2O are shown in 
figures 5-8 for a range of values. The topology of this field (the structure of its foli- 
ation by isosurfaces, or, equivalently, the attractor, repellor, and saddle connections 
of its gradient) appears identical to that of V^p in the case of water, as well as in 
the cases of all other molecules examined by the authors. 

There are some geometrical differences between the two fields. The presence of 
the factor p"^/^ in the expression for Q implies that the magnitude of Q may be large 
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far from the nuclei. In Bohm's picture, this factor accounts for nonlocahty (distance 
independence) in quantum mechanics. In the molecules we examined, this factor 
accounts for structures in Q extending over larger distances than those in V^p. It 
may also be that the quantum force VQ acting on the electrons is responsible, in 
part, for the "concentration" of electronic charge (i.e. regions of Bader Laplacian 
< 0) in bonding and lone pair locations. 

While other authors have compared the structure of V^p with V^i/p/^yp for 
atoms[0| and molecules|]Il[, this is the first time it has been done in three dimen- 
sions, and to our knowledge the first time that a connection of the Bader chemical 
theory to the de Broglie-Bohm quantum potential has been made. 
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figure 7 




figure 4 



figure 8 



figures 1-4: Each figure sliows four views of an isoscalar 
surface of tlic Laplacian of tlic charge density for H2O. 
Figure 1 (V^p <^ 0) shows local minima at the nuclei 
and at the locations of the oxygen lone pairs. Figure 2 
(V^p < 0) shows additional local minima on the bonds. 
In figure 3 (V^p = 0) the outline of the entire molecule 
is visible. Figure 4 (V^p > 0) shows regions of valence 
charge depletion. 



Figures 5-8: Each figure shows four views of an isoscalar 
surface of the (one electron) Bohm quantum potential 
for H2O. The topology and symmetry are identical to 
the Bader Laplacian (V^p, figs. 1-4) for H2O and other 
molecules. 
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